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SUMMARY 



A method -which may he used at high supersonic Mach numbers is 
described for calculating the flutter speed of wings having camber in 
their deflection modes. The normal coupled vibration modes of the wing 
are used to derive the equations of motion. Chord deflections of the 
vibration modes are approximated by polynomials. The wing may have a 
control surface and may carry external stores although no aerodynamic 
forces on the stores are presented. The aerodynamic forces that are 
assumed to be acting on the wing are obtained from piston theory and 
also from a quasi-steady form of a theory for two-dimensional steady 
flow. Airfoil shape and thickness effects are taken account of in the 
analysis . 

The method is used to calculate the flutter speed of seme wings 
which had been previously tested at Mach numbers of 1-5 to 3 - 0 . Com- 
parison of the calculations and experiment is made for flat-plate 60 ° 
and 14 - 5 ° delta wings and also for an un tapered 45° sweptback wing. 


INTRODUCTION 


Tie requirements for low drag at high speed have led to the use of 
thin low-aspect -ratio wings. Structural weights must also be minimized 
to avoid penalizing performance . When such wings are designed primarily 
for the strength required to carry a given static load, the relative 
stiffness is considerably reduced, particularly in the resistance to cam- 
ber deformation. This tendency for wings to deform In camber has intro- 
duced a new element in the flutter picture, namely, a need for a method 
of analysis that can take into account such deformations. Such a pro- 
cedure should be suitable for programing on digital computing equipment 
and could be used as a reference calculation for the correlation of data. 
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In reference 1, several aerodynamic theories were discussed and used 
in comparative flutter calculations. Piston theory, described in detail 
in references 2 and 3, was one of the simpler methods discussed. It was 
shown to give results in agreement with more exact theories when used in 
. two-degree-of -freedom flutter calculations for Mach numbers greater than 
about 2. It is easily adaptable to the present problem of including cam- 
ber deformations, because of its inherent simplicity. A quasi-steady 
aerodynamic method, deduced from second-order steady-flow theory, was 
also tried for comparison with piston theory. Both of these methods 
all ow the effects of airfoil shape and thickness, which are not treated 
in the linearized aerodynamic theories, to be included. 

In this report a procedure for flutter analysis utilizing the nor- 
mal (coupled) vibration modes, which include cambering deflections, as 
the structural input is presented. The first part of the report contains 
the derivation of the flutter determinant. The second part has the 
results of several applications to low-aspect-ratio wings and presents a 
comparison with same experimental results. In order to illustrate the pro- 
cedure, a sample calculation is given in the appendix. 


SYMBOLS 


A iJ ,B iJ ,G iJ' D IJ Burface integrals of wing properties 
a speed of sound, ft/sec 



b 

F 

f ± (x’,y') 


coefficients of polynomials approximating mode shape 
<1 at station n in mode i 

local semi chord, ft 

dissipation function (eq. 11) 

modal function for mode i 


fj-’w 
0 = 


polynomial which approximates mode i 


a _ (y + Dm 4 - 4p 2 
2 p3 


at station n 
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k 

k(x«,y') 

l 


M 




m(x',y') 

“b 



n 

N 


P,Q 

P 


a 

v' v 


T 

t 


structural damping coefficient 
number of modes in analysis 

mass mome nt of inertia of store in pitch, about attach- 
ment point, lb-sec2-ft 

number of polynomials used to approximate chordwise 
mass distribution 

reduced frequency, boi/V 

stiffness influence function, lb/ft/sq ft 

wing semispan, ft 

x- coordinate at ends of mass interval 
Mach number 

generalized mass for mode i 

mass per unit area, slugs/sq ft 
mass of store, lb-sec^/ft 

coefficients of polynomials approximating mass distri- 
bution at station n in interval k 

index indicating spanwise station 

number of spanwise stations used to approximate wing 
properties 

degree of polynomials 
pressure, lb/sq ft 

static unb alan ce of store in pitch about attachment 
point, lb-sec^ 

thickness terms defined by equation ( 26 ) 

kinetic energy 
t ime , sec 
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potential energy 
velocity, ft/sec 


downwash velocity; positive measured away from surface, 
ft/sec 


■ 


+ ig) 


x',y',z' 


Cartesian coordinates 


i< n )(x) 


coordinates along chord, in fraction of local chord, 
and along span, in fraction of semispan, respectively 

coordinate of flap leading edge, fraction of local chord 

function describing airfoil surface, positive away from 
mean surface 

function describing airfoil contour at station n, 
positive away from airfoil mean line 

vertical coordinate of deflected wing mean surface, 
positive up 


ratio of specific heats 

vertical, displacement of deflected wing mean surface at 
store attachment point in mode i 

slope of deflected wing mean surface, relative to 
undeflected mean surface, in stream direction at 
store attachment point in mode i 


b 2 dy‘ 


S ± (t) 


generalized coordinate of motion in mode i 


amplitude of generalized coordinate for mode i 
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p 

CD 

«L 


air density, slugs /cu ft 
circular frequency, radians/sec 
ith natural circular frequency 


Subscripts : 
00 

p ^ q .^,6 

R 

K 


indices indicating mode 
free stream 

indices pertaining to terms in polynomials 
reference 

index -which, denotes the interval along chord over which 
mass distribution is approximated by a parabola 


Superscripts: 

1, 2, 3 } • • • n spanwise station 
Matrix notation: 


[] 

n 


square 

diagonal 

A dot over a symbol indicates differentiation with respect to time. 


ANALYSIS 


In order to develop the procedure, certain classical methods of 
flutter analysis are to be applied. Consider a wing of arbitrary plan 
form whose mass, stiffness, and geometric properties are known. A right- 
hand coordinate system x',y',z', as shown in figure 1, with the origin 
at the leading edge of the wing root is used to describe these properties. 
The wing has a flap which is hinged about its own leading edge. Qhe air- 
foil shape is described by a surface Z(x r ,y r ) which is measured from 
the mean surface of the wing and is taken as positive away from the mea n 
surface. Ihe mass distribution is m(x',y' ), the mass per unit area at 
any point on the wing. Tb* stiffness influence function is k(x',y'), 
an equivalent spring constant at any point on the surface. 
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Derivation of Equations of Motion 

Assume that the deflections can he described by superposing the 
first I vibration modes of the wing. It is assumed for convenience that 
these modes sure normal modes. Then, if fj_(x',y') is the model function 

describing the deflected position of the mean surface in the ith mode and 
^(t) is the corresponding time -dependent generalized coordinate , the 

complete description of the position of the wing is 

I 

z(x',y',t) = ^ (1) 

i=l 

With the position of each point so described, the equations of motion 
can be written with the aid of Lagrange 1 s equation in the form 



where T and U sire the kinetic and potential energies, respectively, 

F is a dissipation function, and is a generalized force. 

Kinetic and potential energies .- The kinetic and potential energies 
of the system can be determined by using the properties of normal modes. 
The kinetic energy will be 

T = i JJ m ( x '»y , )|*(x , ,y , ,t)J s*' d y' (3) 

s 

where S denotes integration over the entire surface of the wing. Com- 
bining equation (l) with equation (3) yields 

11 

1 - i I Y. iT m( xSy'jffiCxSy 1 )] [fj(x',y')]dx' dy' 

1=1 J=1 S 

00 

Since the modal functions are normal modes, the orthogonality relation, 
derived, for example, in reference 4, applies . That Is, 
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<4 


jJ'm(x',y)[f 1 (x' J y)][fj( x',y' )] dx' dy' = 

S 

JJ'm(xSy')|f 1 (x>,y’)] [fj(x*,y')]dx' dy- = 0 


S 


(i = 3) 
(i t 3) 


(5) 


where Is the generalized mass In mode I. Thus, the kinetic energy 

can he written hy using equation (5) as 


T = 



1-1 


( 6 ) 


Die potential energy of the system Is 

U = i JJ k(x' ,y 1 ) [z(x* ,y ' ,t)J dx- dy' (7) 

S 

Substituting equation' (l) into equation (7) yields . 

I I 

-■III s i (t)s j jfj(x',y')Jdx' dy- (8) 

i=l 3=1 S 

Now, the orthogonality relation of equation (5) can he written as 


JJ k(x' ,y ' ) Jf^Cx' ,y' )J |f ^ (x- ,y- )j dx' dy- = (l = J) 

8 

JT Hx * ,y« ) [f i(x- ,y- )] [f J (x- ,y - )] dx' dy- = 0 (i / J) 

S 


(9) 


where 
energy Is 


is the natural frequency of mode i. Thus, the potential 


U 








( 10 ) 
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Dissipation function .- Reference 5 discusses the development of a 
dissipation function which is convenient for allowing for internal damping 
in flutter analyses for the particular case of harmonic time variation. 
This function is a convenient expression which, when differentiated with 
respect to the velocities ^(t), gives the generalized damping force in 

the system. For a wing which has structural damping, the dissipation 
function can he written 



i=l 


(H) 


where g^ is the structural damping coefficient in mode i. 

Equations of motion .- Operating on equations (6), (10), and (11), 
as indicated hy equation (2), yields the equations of motion 

MjV^t) + g^ i ± (t) + ox^^t) = £ Qj (12) 

J=1 


Generalized Forces 

The external forces acting on the wing which influence the flutter 
system are aerodynamic forces. These aerodynamic forces will he calcu- 
lated hy two procedures: (a) piston theory and (h) a quasi-steady method 

based on second-order two-dimensional theory. 

Piston theory .- Piston theory has heen developed and discussed in 
references 1 to 3. The basic assumptions of this theory are slender 
profiles and high Mach numbers (m 2 » l.o). As a result of these assump- 
tions, the x' component of fluid velocity changes very little along the 
profile and is always much greater than the speed of sound. This means 
that disturbances at one point on the wing induce only a s mall effect at 
another point. Piston theory neglects these sma ll induced effects by__ 
assuming a point.. function relation between pressure and downwash velocity. 
This relation assumes that the local pressure on an airfoil is related to 
the local fluid velocity normal to the free stream in the same m a n ner 
that the pressure on a piston in a one-dimensional channel is related 
to the velocity of the piston. Obviously, three- dim ensional effects are 
not included in piston theory. However, these effects on wings are 
relatively smal 1 at high flight speeds . 

The second-order piston-theory relation between pressure and down- 
wash, as given in reference 3, is 
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Die downwash (or piston velocity) at the point under consideration is w 
and w/a should he less than 1.0. Hiis expression has been found to 
apply reasonably well for Mach numbers greater than about 2.0. 


Quasi-steady second-order theory .- An expression similar to piston 
theory for the pressure on a wing can be obtained from the second-order 
theory for steady flow which is given, for instance, in reference 6. If 
the local steady angle of attack at a point on the surface is replaced 
by the instantaneous (slowly varying, hence quasi-steady) local angle of 
attack, the pressure at the point is 


P - P 


00 



w 

a 


+ 



+ 1) - ^P 2 /wN 2 




Equation (l4) approaches the piston- theory pressure (eq. (13)) if the 
Mach number Is very large; that is, £ *** M. Equation (14) contains no 
as s ump tion about the Mach number being large; this expression could give 
more useful results for pressure than piston theory in a range of Mach 
numbers around 2.0. 


Downwash .- Hie downwash appearing in equations (13) and (l^) Is the 
Instantaneous local slope in the stream direction multiplied by the free- 
stream velocity. It can be separated into displacement and airfoil- shape 
terms (as in ref. 2) which ere not the same for the upper and lower sur- 
faces of the wing. If z(x , ,y , ,t) is the position of the mean surface 
and Z(x' ,y* ) describes the contour of the airfoil surface measured from 
the mean surface, the downwash will be 


For the upper surface: 


v - (v + |-')z(x , ,y , ,t) + V > v I (15a) 

\ ox' ot/ dx* 

For the lower surface: 

v --( v a? + to) 

Downwash is positive when measured away from the surface of the wing. 

Differential pressures Pressures acting on the wing surface are 
found by substituting equation (15) into equation (13) or (l4). The 
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differential pressure, in the positive z' direction, acting across 
the wing is then determined by subtracting the pressure on the upper 
surface from the pressure on the lower surface. If only time -dependent *■ 

terms and terms which are lin ear in displacement are retained, the dif- 
ferential. pressure when the piston theory is used is 


Ap(x',y',t) = -2pajl + G Z(x',y') ^- + |^jz(x' ,y' ,t)j (l6) 

*y 4 _ ~| 

where G = M — . Use of the quasi-steady theory would multiply the 

right-hand side of equation (l6) by a factor of M /jj and change the fac- 
- _ M^(y + 1) - kfi 2 


tor G to G = 
large Mach numbers 


2P' 


The factor G approaches G for very 


The modal representation of the wing's motion (eq. (l)) can be sub- 
stituted into equation (l 6) and the resulting expression for the differ- 
ential pressure is 


I 

Ap(x',y',t) = -2pa^ t + G 
J=1 L 


7 Z(x*,y*)jjv ) ij(t)+fj(x',y')ij(t)j 


(IT) 


Generalized aerodynamic forces .- Ihe generalized aerody nam ic force 
in mode i cairbe interpreted in terms of the total virtual work in 
mode i due to the aerodynamic forces. As shown in reference k, this 
generalized force is determined by the integration over the surface of 
the product of the differential pressure distribution and the deflection 
distribution in mode i. The result for piston theory, when equation (IT) 
is used for the pressure, is 




j=l s 


fj(x',y')ij(t) 


) ^(xSy'jJ 


dx' dy' 


( 18 ) 


For systematic evaluation of equation (lS), it is convenient to separate 
the terms as follows: 


l- 
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^ ' _2pa I [ V ( A U + Q °l3) S J (t) + ( B H + ®ij)«3 {t >] (19 > 


J=1 


where 




ppVbt .(x^y*)"],- 

■ 11 ax- ] k ( *'’ y 'i] to ' ^ 

s 


B. 


U 


'U 


D<* = 




M jUSy')] ^(x'^y* jja*' ay' 

1 |sF z <*' > y ' >] |f i<*’ - y ' >] ta ' *' 

s 

JJ j^T Z(x« ,y' )J [f J (*• J,y * )] [f^x* ,y' )J dx« dy’ 


s 


( 20 ) 


Equations (20) are parameters ■which are constant for a given wing, 
depending only on the mode shape and. the contour of the wing surface. 
The quasi-steady aerodynamic theory will alter the generalized force 
expression (eq.. (19)) "by the multiplicative factor M/p and substitute 
G for G. 


Flutter Determinant 

Bae flutter determinant for this system may be obtained from the 
equations of motion (eq.. (12)) and the generalized forces (eq.. ( 19 )). 

In order to obtain a flutter solution, the assumption is made that the 
damping coefficient in all modes is the same; that is, g-, = gg = • • • = g. 

The assumption of simple harmonic motion is not required when the simpli- 
fied aerodynamic theories are used since the roots of the characteristic 
equation could be examined directly. However, this procedure is presented 
in the usual form for flutter calculations by letting the time variation 
be a simple harmonic variation to obtain the borderline flutter condi- 
tion. If equations (12) and (19) are combined and the terms rearranged, 
the equations of motion at flutter for piston-theory aerodynamics become 
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L J J=lL 


M(’ 


j + GD ± j j j Sj = 0 (subscript i = 1, 2, ... I) (21) 


where 


S(t) 



X = (S) 2 (1 + lg) 

" *R $ 

The subscript R refers to any convenient reference station. 


( 22 ) 


The requirement for a nontrivial solution to equations (21) is that 
the determinant of the coefficients of l ± must vanish. Thus, in matrix 

notation. 



where 



indicates a square matrix and 



indicates a diagonal 


matrix. This equation must be satisfied at any reduced frequency, Mach 
number, and density as long as the motion is simple harmonic. A common 
method of finding the borderline case of neutral stability is by solving 
for the velocity at which zero damping is required to maintain harmonic 
motion . 


This flutter dete rmina nt has been derived by using piston-theory 
aerodynamics. If the quasi-steady pressures given by equation (14) had 
been used, the flutter determinant would be the same except_that 0 
would replace M and the factor G would be replaced by G. 

Ihe flutter equations as derived herein contain airfoil-shape terms. 
The usual assumption of linearized supersonic theory is that these terms 
are negligible . In order to drop thickness terms from this calculation 
and thus have a high Mach number low-frequency approximation for linear 
theory, it is only necessary to set the factor G (or 5 ) equal to zero 
wherever it appears in the flutter determinant-. - 
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Representation of Wing Deflections 

Closed-form evaluation of the surface integrals defined by equa- 
tions (5) and (20) Is not practical for the complex modal deformation 
shapes normally possessed by wings. In this report, a numerical evalua- 
tion of these Integrals is accomplished by approximating the wing charac- 
teristics at a limi ted number of spanwise stations, performing the chord- 
wise Integrations, and then summing spanwise. In order to use this 
method, the wing is divided into N spanwise stations, parallel to the 
airstream. At each of these stations, the chordwise distribution of dis- 
placement in mode i is approximated by a polynomial of degree P over 
the chord ahead of the flap and by a polynomial of degree Q over the 
flap, as shown in figure 1. That is, the total deflection at station n 
in mode i is given by 

A ■ E 

p=0 

■ i (".‘-V 

q.=0 



where x is a nondimens ional streamwise coordinate havi n g the value 
zero at the leading edge and 1.0 at the trailing edge. The coeffi- 
cients (^i^) are conB ' fcan ' fcB - This series of equa- 

tlons approximates the modal function f^(x l ,y ! ) where the y 1 varia- 
tion has been replaced by an approximation at N spanwise stations, 
denoted by the index n. 

The surface integrals of equation (20) can be written in terms of 
the deflection polynomials by substituting equation (24). The integrals 
of equation (20) involve the products of two polynomials or the product 
of a polynomial and the derivative of another polynomi a l. The spanwise 
coordinate becomes y, which is zero at the wing root and unity at the 
tip. The chordwise integrations for Cj_j and Dj_j are done by parts 

in order to get the airfoil cross-sectional properties in convenient 
form. Thus, equations (20) become 
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*« • •/„ W Z •W'-’lH" * cl 1 4/10; 

tj- 0 rrf J- tf^3 


x^laxl ay 


z(-/" i )i i iH“* c K ^(“/loH^V 

p-0 r-0 1 q.-0 «^> 


«j «i 


* ■ */: 1 1 W1W * 1 1 4/l(0.(£4 * 


lO a-0 


*-»*•/: 4 1 zf/io,(^)*? fo.owV 


p-0 3M) 


1-0 r-0 


The t v ^ and s v ^ n ) terms appearing in these equations are constants 

which are determined by the airfoil thickness and cross section and are 
defined by the following equations: 

. (n) / Cz (n> (*i)l .. r x 1 _«tw W l 1 


(n) , N v[z (n, (*l)l r * 1 v-lfz (ll> C*)l. 

'WL-^J- v j 0 * N“ J 


where v takes on the values of the indices in equation (25). The 

quantity Z^(x), illustrated in figure 2, describes the airfoil contour 
at station n and is taken as positive when measured away from the mean 
chord. 

When the thickness goes to zero, the t ^ and s v ^ n ^ terms 

become zero. However, as previously mentioned, a zero-thickness calcula- 
tion can be accomplished numerically by setting G or G equal to zero. 
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Generalized Mass 

Distributed mass .- The generalized mass (defined by eq.. (5)) 

can be computed In any convenient manner . However , If the wing is of 
solid construction or has a»mass distribution which is not too irregular, 
the mass distribution can be approximated at Idle N stations discussed 
in a previous section by a series of parabolas; thus the generalized 
mass can easily be evaluated. If the wing chord at a spanwise station n 
is broken into K arbitrary intervals, as shown in figure 3 , the mass 
per unit area in any one of the intervals k is given by 

>> (X) = (^) K * (»x (n) 1* * (Vi *«**.) 

If the wing has an aileron, the appropriate mass interval can be made to 
coincide with the aileron. That is, will equal x^ for the 

example illustrated . With the mass distribution approximated in this 
manner and the deflections approximated by equation (24), the generalized 
mass from equation (5) becomes 



External stores and pods .- In order to account for rigidly attached 
external stores or pods, additions must be made to the previous ana lysis . 
For the present case, the assumption is made that aerodynamic forces 
acting on the body are negligible. Pit s means that the terms j , 

, and Dj_j defined by equation (20) are unaffected. The only way 

that the stores can enter the problem is through the kinetic and potential 
energies and the dissipation function. 
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If the store is represented by a concentrated mass, inertia, and 
unbalance about its attachment point, the attachment point displacement 
and slope in the stream direction for mode i being given by 5^ 

and respectively, the increment in kinetic energy which must be 

added to equation (4) is 

11 r 

^ = § X X 'i (t)i j (t) kvj + \(\ & 3 * Vi) + V^l (29) 

i*l j=i L -» 

The orthogonality relation between normal modes will apply for modes 
determined with the additional mass of the store in the system. Thus, 
the generalized mass M^ for a wing carrying external stores would be 

Mi - JJ mCx^y'J^CxSy')! dx' dy» + + ^cts 6 ! 6 ! + ^ b « 1 ^|(30) 

S 

With this new definition of the generalized mass, the kinetic and poten- 
tial energies and the dissipation function, appear as before in equa- 
tions (6), (10), and (ll). Therefore, if a wing is carrying stores, 
the generalized mass given by equation ( 28 ) must be increased by an 
increment 


AM ± = m s8l 2 + SB a S ± * t * I^ 2 (31) 

This is the mass effect of the store. If the aerodynamic forces are 
known, they can be included along with the other aerodynamic terms. 


APPLICATION AND COMPARISON WITS EXPERIMENT 


The method described in the previous section has been used to cal- 
culate the flutter characteristics of several cantilevered wings at 
supersonic speeds. Atypical results of these calculations and compari- 
sons with experiment will be presented for 'three wings. Piston- theory 
calculations for Mach numbers as low as 1.3 were made for comparison 
with exper ime nt and with quasi-steady theory, although the piston-theory 
region of applicability is usually considered to be Mach numbers greater 
than about 2. 
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Two of t3ie wings covered by these calculations were 60° and 45° 
flat-plate delta wings . Their airfoil sections were flat with a slight 
bevel at the leading edge and were approximately 0.5 percent thick at 
the root. A third wing studied was a uniform untapered 45° swept wing 
with a panel aspect ratio of 1.33 ■ Its airfoil section was flat sided 
with both edges beveled and was 1.4 percent thick. The geometric, mass, 
and vibration characteristics of these wings are given in reference 7* 
Their experimental flutter characteristics were reported in reference 8. 


Approximation of Wing Properties 

For the wings under consideration, the experimentally determined 
camber mode shapes from reference 7 were approximated by linear, quad- 
ratic, and cubic polynomials. Thus, the maximum value of P from 
equation (24) was 3* None of the wings had flaps; thus x-^ in equa- 
tion (l) may be set equal to 1.0. The accuracy of the approximation 
for camber mode shapes is illustrated in figure 4 where the camber modes 
and their polynomial approximations are plotted at various span stations 
for the first three vibration modes of the 43° delta wing. Figure 4 
also shows the large amount of camber deformation present in wings of 
this type. The approximations are seen to be best near the tip where 
the deflections are the greatest. Also, the approximations are better 
in the lower modes (figs. 4(a) and 4(b)) than in the third mode 
(fig. 4(c)). Better approximations in the higher modes could probably 
be obtained by using polynomials of higher degree. 

For this group of wings, ten spanwise stations (N * 10) at 95- * 

85-, 75-, • . . 5-percent span were chosen for computing the system 
properties. Spanwise integrations were made by using a ten-point 
numerical integration scheme. 

The 60° and 45° flat-plate delta wings had constant mass distribu- 
tions, the bevel being neglected. The 45° swept wing was a plate with 
beveled edges such that the mass distribution was proportional to the 
thickness at any point. The chordwise mass distribution was constant 
over the center portion of the chord and, had a linear variation over 
the beveled edges. Equation (27) described this distribution accurately 
with K = 3- 


Results of Flutter Calculations on Flat-Plate Delta Wings 

The calculations on one wing, the 60° flat-plate delta, are fully 
described in the appendix. The results are shown in figure 5 and com- 
pared with experiment from reference 8. Four separate results of this 
three-mode calculation are shown: for piston theory with and without 

the airfoil- thi c kne ss term and for quasi-steady theory with and without 
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the thickness term. The flutter frequencies, shown as ratios to the 
lowest natural frequency, are calculated by all methods to within 15 per- 
cent for Mach numbers from 1.3 to 3-0. The flutter boundary is shown in 

figure 5 as a plot of the stiffness-altitude parameter against 

Mach number. 

The stiffness -altitude parameter, used in figure 5 and subsequent 
figures, depends only upon the physical properties of the wing (in 
particular, the stiffness) and upon the atmosphere in which it operates. 
Its value increases as either the altitude or the stiffness increases. 
When plotted as the ordinate against Mach number, as in these figures, 
constant dynamic pressure curves will appear as radial lines through the 
origin. The stable region will be above the flutter boundary. 

In this case, piston theory predicts flutter on the conservative 
side within 23 percent of experiment at M. = 1.3 and within 9 percent 
at M » 3*0. Quasi -steady theory predicts flutter still more conserva- 
tively, within 45 percent at M => 1.3 and within 12 percent at M = 3.0. 
In both instances, as would be expected for a thin, flat-plate wing, 
thickness terms in the aerodynamic theory had very little effect. It 
should be noted that, on a percentage basis, both theories give better 
agreement with experiment as the Mach number becomes larger and the 
reduced frequency becomes smaller. 

The results of similar three-mode calculations on a 45° flat-plate 
delta wing are shown in figure 6. The experimental data are from refer- 
ence 8. Agreement between calculation and experiment is better than that 
for the 6o° delta wing but the calculations are now slightly unconserva- 
tive. Piston theory deviated from experiment by 14 percent at M = 1.3 
and by 3 percent at M = 3} thus, the accuracy improved as the Mach num- 
ber increased. Quasi-steady theory showed better agreement than piBton 
theory in this case, 2 percent at M = 1.3 and M = 3* Again, thick- 
ness terms in the aerodynamic theory had very little effect on the cal- 
culated flutter speed. Calculated flutter frequency was good, within 
6 percent at all Mach numbers . 


Results of Flutter Calc ula tion a on a 
Flat-Plate 45° Swept Wing 

The results of calculations made with three modes on the untapered 
45 swept wing are shown in figure J. Agreement with exper ime ntal data 
from reference 8 is not nearly as good as that for the delta wings. The 
calculated flutter frequency is seen to be high by about 95 percent for 
all four methods of calculation. The flutter bo undar y is predicted uncon- 
servatively by all methods. A possible explanation for the sizable dis- 
crepancy between calculations and exper ime nt is the large tip effect. 
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Especially at lower Mach numbers, this wing will have a large portion of 
its surface influenced by the tip, which is parallel to the flow. Hiis 
effect was not present on the delta wings. 


CONCLUDING REMARKS 


A high supersonic Mach number flutter-calculation procedure for 
wings with camber in their vibration modes has been described. The nor- 
mal (coupled) vibration modes of the wing are utilized and their camber 
deflections are approximated by polynomials. The wing may have a con- 
trol surface and may carry external stores. Aerodynamic forces on the 
stores are neglected but could, be included if desired. The aerodynamic 
forces acting on the wing are obtained from piston theory and also by a 
quasi-steady method based on second-order supersonic theory. Both aero- 
dynamic theories allow airfoil shape and thickness terms to be included. 
Use of the piston theory requires the Mach number squared to be much 
greater than 1 and both aerodynamic theories require the reduced fre- 
quency to be low. 

The method was used to calculate the flutter speed of three canti- 
levered wings from Mach numbers of 1.3 to 3-0. Good agreement with 
experiment was found for a 6o° and a 45° flat-plate delta wing. Sizable 
differences from experiment resulted when the method was applied to an 
untapered 45° swept wing. 


Langley Aeronautical Laboratory, 

National. Advisory Canmittee for Aeronautics, 
Iangley Field, Va., August 8, 1958- 
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APPENDIX 

** 

SAMPIE CALCUIATIONS ON A 60° FIAT-PLATE DELIA WING 


The method of flutter calculations described in the body of the 
paper Is illustrated here for the 60° flat-plate delta wing. The wing 
was made of a flat plate of metal and was cantilever mounted at the root. 
All computations were performed on digital computers. 

The modal data needed in the calculation are given In table I. The 
first three vibration modes are used and the wing deflections are deter- 
mined at' ten spanwise stations. In order to specify the chordwise poly- 
nomial for each vibration mode and at each station n, the degree of_ the 
polynomial P and the coordinates of the polynomial at P + 1 chordwise 
points are given. Then, by solving a set of P + 1 s imult aneous equa- 
tions, the constants of equation (24) can be evaluated for each sta- 
tion n and each vibration mode. 

With the constants of the deflection polynomials determined, the 
flutter-determinant elements given by equations (25) and (28) can be 
evaluated by using the mass and geometric data given in table I. Note 
that for this flat-plate wing the mass distribution is constant. The 
airfoil has been assumed to have zero thickness at the leading edge, 
flat sides, and nonzero thickness at the trailing edge; thus all thick- 

(*i) 

ness terms are zero except tg • The wing does not have a flap so 
the dimension is set equal to 1.0. The elements for the flutter 

determinant computed by using this data are 


0 0 

2.2970845 0 

o 0.92344991 


Kl 


10 


-4 


1.1087469 

0 

0 
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b E 2 p 4 .jj ■ 10 2 

0.24745638 
-o . 24950129 
0.068579763 

1.7161381 

-0.15541714 

- 0.24973402 

-0. 47849789" 
0.34179249 
0.20163600 

[ b ij] - 10-3 

” 3.9726354 

-0.24698926 

-0.14509758 

- 0.24698912 

8.2304250 

0.071423539 

-0.14509758 ~ 
0.071423468 
3.3087173 

- w 4 

’ 0.50658132 
-0.25845494 
0.42135938 

1.1599243 

1.2900028 

1.1803951 

- 0.014184357 

-0.34968308 

- 0.037120715 

”R [B ±J j = m -5 

~ 3.1259500 
-0.50451053 
_ 2.7703439 

- 0.50451053 
6 . 9119095 
0.54545789 

2.7703439 " 

0.54545789 

2.6712828 


(Al) 


These el em ents, computed from equations (25) and (28), together 
with the natural frequencies, are fixed for this wing. The aerodynamic 
parameters, Mach number and density, can he varied Independently In the 
flutter determinant of equation (23) . For example. If the wing fluttered 
at a density of O.OOO 63 slug/cu ft at Mach number 3 . 0 , the factor G 
Is 3*6 and the flutter determinant Is 
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a.asx - U + o-ooBfijasBrn^j ♦ <h«*onn&^ - ajcMfaooo^lj 

- 11 rrrKHOTno^ 

1 

-o-mWSJUST^j - - o-xe^rn . cumeegrTO^Jj^ * o-uuanga»i^Lj 

CJxgsaw^j * ojxnosoagoa^) 

' i 

- 0 

Mooewin^ji - o-acrtaBeMD^j * o-corngnatmus^ 

%toBi n - ojBgesm * oAmtanb^^ * ojtawws^y 



Iblg Cmyj cX uetcgsliittut Ccujl be Solved fur flutter Sjjccd cuiu. fjre^ulcuCjr 111 a in. Hub ex uf wcvjr d • 

Qhe moat ccjanon method la to substitute successive reduced frequencies and to solve each 

tine for the complex root X, defined by equation (22) . At flutter the Imaginary part of this 
root, proportional to the damping required to produce harmonic motion, la zero. Plots of a amp -trig 
coefficient and frequency obtained In this manner are shown In figure 8. From this figure It Is 

SeeH ULuib 'the TOOt C0rfe6p0mfli5g tu ‘the Second uatUTul Vlumtiun SOde ueCdueS Uu5~tabXs at 

j = 53 and a frequency which la 2-53 times the first natural frequency. Thus, the flutter fre- 
quency is 156 cps and the flutter velocity is 1,842 ft/sec. This result is plotted in figure 5- 

If It were desired to use quasi-steady aerodynamic s at the same Mach number , B = 2.83 
would be substituted for M and G = 3.58 for G in equation (23) • The terms G or 5 can 
be set equal to zero for a zero-thickness calculation. 


1 i 


Mi VGVH 
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(a) First natural mode. 

Figure 4.- Comparison between measured chord deflection mode shapes and 
their polynomial approximations for the first three natural modes of 
a 4-5° delta wing at various span stations. 
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(c) Eiird natural mode. 
Figure 4.- Concluded. 
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